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Abstract 


A  perturbation  analysis  was  done  on  a  satellite  in 
orbit  about  L3.  A  modal  control  scheme  for  linear  systems  was 
then  applied  to  the  linear  and  nonlinear  perturbations. 
Rectification  of  the  satellite  orbit  was  achieved  by  Encke's 
Method.  The  effects  of  the  moon's  inclination  and  eccentri¬ 
city  on  the  orbit  stability  was  investigated.  The  controller 
stabilized  the  linear  and  nonlinear  perturbations,  with  the 
average  control  acceleration  being  of  the  same  order  of 

magnitude  as  that  for  earth-synchronous  satellites. __ 
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MODAL  CONTROL  OF  A  SATELLITE  IN  AN 
UNSTABLE  ORBIT  ABOUT  L3 


I.  Introduction 


Background 

The  need  for  total  global  reconnaisance  via  sensor 
bearing  satellites  has  been  obvious  since  Sputnik  I. 
Equatorial  synchronous  satellites  provide  only  partial 
coverage,  not  being  able  to  "see"  polar  activities. 
Synchronous  satellites  also  require  station  keeping  devices, 
e.g.  reaction  jets,  since  they  are  in  an  inherently  unstable 
orbit  (ref.  1).  The  cost  of  this  control  in  terms  of  weight 
may  be  excessive,  to  the  point  of  jeopardizing  the  basic 
mission  of  the  satellite.  If  a  satellite  can  be  put  into  an 
unstable  orbit  and  kept  there  with  a  minimum  amount  of 
control,  then  less  weight  need  be  allocated  to  the  control 
system.  Therefore,  the  need  of  minimizing  control 
requirements  becomes  clear.  There  also  arises  the  need  for 
a  point  or  points  in  space  in  which  a  satellite  can  be 
placed  and  stabilized  with  a  minimum  of  control  cost  and  be 
able  to  provide  total  global  coverage. 

There  are  points  in  space  where  an  object,  when 
placed,  will  tend  to  remain  if  left  unperturbed.  In  the 
Earth-Moon  system  these  points  are  known  as  Lagrangian 
points  Ll  through  L5  (See  Figure  1).  These  points  arise  in 
the  Earth-Moon  system  for  the  restricted  three-body  problem. 
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The  model  is  such  that  the  Earth  and  Moon  rotate  about 
their  barycenter  in  circular,  planar  orbits.  The  collinear 
points  LI,  L2,  and  L3  are  inherently  unstable,  and  an  object 
placed  at  these  points,  when  perturbed,  will  tend  to  leave 
the  vacinity  of  these  points.  The  points  L4  and  L5  are 
inherently  stable  and  an  object  placed  at  these  points  will 
tend  to  remain  even  if  perturbed.  But  when  a  fourth  body  is 
considered,  the  sun,  all  these  points  become  unstable. 

It  is  known  that  periodic  orbits  about  these 
equilibrium  points  exist  (ref.  2).  In  the  restricted 
three-body  problem,  if  the  ratio  of  the  two  larger  masses  is 
less  the  1/27,  orbits  about  the  triangular  points  can  be 
shown  to  be  unstable  (ref.  3).  Wiesel,  in  an  earlier  work, 
established  an  orbit  about  L3  for  the  restricted  four-body 
problem  (ref.  2).  Figure  2  shows  this  orbit.  This  orbit  is 
found  to  be  unstable.  Therefore,  a  means  must  be  found  of 
stabilizing  a  satellite  placed  in  this  orbit. 

The  Lagrange  points  lie  in  the  Earth-Moon  plane  and 
not  the  equatorial  plane  which  makes  polar  reconnaisance 
from  a  satellite  in  orbit  about  L3,  L4,  and  L5  possible.  It 
is  the  lack  of  total  global  coverage  by  Earth  synchronous 
satellites  that  makes  the  prospect  of  putting  one  satellite 
in  orbit  about  each  of  the  Lagrangian  points  L3,  L4,  and  L5 
desirable.  However,  it  must  be  done  such  that  control  costs 
are  not  excessive. 
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Problem  and  Scope 

In  an  earlier  work,  Shelton  developed  a  modal 
control  scheme  for  time  periodic  linear  systems  (ref.  4). 
This  scheme  was  then  applied  to  a  satellite  in  an  unstable 
orbit  about  the  Earth-Moon  Lagrange  point  L3.  With  the 
controller  operating,  the  orbit  becomes  a  stable  reference 
solution  (trajectory)  and  classical  perturbation  techniques 
could  be  used.  Shelton's  controller  was  designed  based  on  a 
model  which  neglected  the  free  lunar  eccentricity,  lunar 
inclination,  and  the  eccentricity  of  the  sun.  When  these 
effects,  particularly  the  free  lunar  eccentricity,  were 
included  in  the  model,  the  controller  could  no  longer 
stabilize  the  orbit. 

Using  Shelton's  periodic  orbit  as  a  stable  reference 
solution,  Ehrler  (ref.  5)  did  a  first  order  perturbation 
analysis  of  the  effect  of  the  free  lunar  eccentricity.  With 
the  model  changed  to  include  the  Moon's  free  lunar 
eccentricity,  Shelton's  controller  responded  to  both  the 
original  instability,  which  it  had  previously  stabilized, 
and  the  purely  oscillatory  linearized  perturbation  from  the 
free  lunar  eccentricity,  which  in  combination  it  could  not 
stablize.  This  led  to  Ehrler  splitting  the  perturbation 
solution  into  its  free  (homogeneous)  and  forced  (particular) 
solutions  which  is  the  basic  theory  of  linear  differential 
equations  with  a  forcing  function.  Ehrler  then  modified  the 
controller  to  suppress  only  the  "free"  components  causing 
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the  original  instability.  However,  when  the  Moon's 
eccentricity  was  included  in  the  system  model,  the 
perturbation  was  so  large  it  drove  the  satellite  outside  the 
linearized  region  the  controller  was  designed  for,  rendering 
it  unstable. 

Instead  of  linearizing  the  perturbation  equation,  we 
can  choose  to  redefine  the  perturbation  so  that  it  is  split 
into  two  parts,  the  first  consisting  of  the  original 
instability  (linearized  and  controllable),  and  the  second 
consisting  of  the  perturbation  from  the  free  linear 
eccentricity  and  the  nonlinear  terms  of  the  original 
instability.  By  subtracting  off  the  second  part  from  the 
total  perturbation  we  can  modify  the  controller  to  suppress 
only  the  free  component.  This  will  hopefully  be  successful 
in  stablizing  the  satellite  in  the  reference  periodic  orbit, 
even  if  the  perturbation  is  so  large  as  to  drive  it  outside 
of  the  linear  range  of  the  controller  design.  This  is  the 
goal  of  this  study. 
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II.  Problem  Analysis 

Coordinate  System 

The  coordinate  system  used  is  shown  in  Figure  3. 

The  coordinate  system  rotates  with  the  same  angular  velocity 
of  the  moon.  The  origin  of  the  system  is  located  at  the 
center  of  mass  of  the  Earth-Moon  system,  the  barycenter. 


Dynamics 

In  an  earlier  work,  wiesel  developed  the  dynamics 
for  the  restricted  four-body  problem  (ref.  2).  The 
Hamiltonian  was  found  to  be 


H  - 


i(p2  +  p2  +  p2)  +  (p  R  -  p  R  ) 
1  Fx  *y  *x  y  *y  x 


m. 


M 


p  ( r  *  -A)  )  +  p  (r  •  +A)  ) 

px  u'py  *y  py  *^p* 


Gm. 


Gm, 


Gm„ 


(R-p) 


m„ 


r+R 


ml+m2 


m. 


R  - 


ml+n*2 


(1) 


The  standard  units  and  conventions  of  the  restricted 
three-body  problem  have  been  adopted  here  (See  Appendix  A). 


Equations  of  Motion  and  Equations  of  Variation 
The  state  vector  formulated  is 


X  = 


i  -  1,  2,  ...  n 


(2) 
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Where  the  q^'s  represent  the  generalized  coordinate  and  the 
p/s  the  associated  conjugate  momenta  in  the  rotating  Earth- 
Moon  frame.  The  equations  of  motion  are  found  from  Hamilton's 
equation 


*i  - 

a  h 

Pi  = 

-  2  H 

(3) 

a  Pi 

This  can  be  put  into  the  form  x  = 

[Z]  *H 

(4) 

where  Z  is 

a  skew-symmetric 

2nx2n 

matrix,  n 

being  the  number 

of  degrees 

of  freedom.  Z  is 

of  the  form; 

[2]  = 

f  to]  [In ] ~] 

(5) 

[-In]  [o] 


Here  [In]  is  the  nxn  unit  matrix. 

In  a  similar  manner  the  equations  of  variation  can  be  found 
from 

<&  =  [z]  fH]<&  (6) 

Where  ['fr]  is  the  2nx2n  Hessian  matrix  whose  elements  are 
“Hij  =  32H  i,j  -  1,  2n 

(ref.  6:230). 

We  can  put  these  time  periodic  variational  equations 
in  the  form 

ix  =  A(t)  £x  (7) 

The  periodicity  is  contained  in  the  coeffiecient  matrix  A(t) 
such  that 

A(t)  *  A(t  +  T)  where  T  is  the  period  of  the  system. 
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III.  Floquet  Solution  For  Periodic  Systems 

For  a  time  periodic  system  such  as 

£x  *  A(t)ix  (8) 

There  is  a  fundamental  matrix  or  state  transition  matrix 
$i(t,  t  ),  such  that 

&x(t)  =  $(t,  0)  £  x(0)  (9) 

The  fundamental  matrix  is  a  function  of  only  the  initial  and 
final  times  and  has  the  property 

|<t,  tQ)  =  $(t,  tx)  J(t1#t0) 

The  n  columns  of  $(tft0)  each  satisfy  the  differential 
equation 

4>(t)  =  A(t)^>(t)  (10) 

If  §(tQ,  tQ)  is  the  identity  matrix,  i.e. 

i'V  to»  ■  1 

Then  §(t,  tQ)  is  the  principal  fundamental  matrix. 

Floquet  Theory  states  that 

§(t,  T)  =  <§(t)C  (11) 

where  C  is  the  monodromy  matrix.  Substituting  t  *  0  into 
Equation  (11)  we  see  that  §(T)  ■  C  (12) 
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By  investigating  §(T)  we  obtain  the  eigenvalues  and 
eigenvectors  at  t  >  0.  He  can  obtain  §(T)  by  using  appro¬ 
priate  initial  conditions  and  integrating  Equation  (10)  for 
n  columns  of  §,  over  the  period  T.  Ploquet's  Theory 
states  that  the  general  solution  of  Equation  (8)  is  of  the  form 


k 


x(t) 


n 

Y.  ‘i  e 

j=i 


(13) 


where  8^  and  are  constants  and  the  Ej  are  periodic 
functions  with  periodic  (T).  The  n  X..  can  be 
arranged  in  a  matrix  of  Jordan  canonical  form,  i.e. 


J  = 


\ 


\ 


\ 


\ 


'n 


(14) 


(ref.  9). 


The  Xj  are  called  the  characteristic  exponent  or  Poincare 
3  \,t 

Exponent  and  the  quantity  e  J  termed  a  characteristic 
multiplier  of  A(t). 

In  order  to  compute  E^  and  \j  Shelton  assumed 
that  he  needed  to  excite  one  mode  of  the  system  described  by 


(8).  Letting  Bj  * 
from  (13),  we  get. 


0  for  j  ^  i  and  8^  =  1  for  i  =  j 


(—  -  Xft 

ox(t)  *  E^t)  e 


(15) 
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Equating  (9)  and  (15)  we  obtain 

Ei(t)  e  =  §(t,  0 ) ox ( 0 )  (16) 

=  $(t,  0)  Ea(0)  (17) 

E-^t)  is  periodic 

Hence  E(T)  =  E^(0),  which  results  in 
1-  -  XiT 

ox(T)  =  Ei(T)  e  (18) 

X-T 

=  Ei(0)  e  A  =  ${t,  0)  E^O)  (19) 

This  leads  to  the  eigenvalue  problem 

T  Xi<T) 

($(T,  0)  -  (e  )I}  Ea(0)  ■  0  (20) 

Here  we  see  that  the  monodromy  matrix  has  eigenvalues  eXT 
and  its  eigenvectors  are  E^(0).  To  obtain  the  complete 
solution,  E^(t)  must  be  computed  over  the  period. 

Determining  E^(t)  over  one  period  determines  E^(t) 

for  all  time.  Letting  t  =  0,  Floquet's  Theorem  also  states 

that  we  can  describe  the  fundamental  matrix,  §(t,  0),  as 

|<t,o)  =  _/\_(t)eJt_/\__1(0)  (21) 

where  _f\_  is  the  matrix  of  eigenvectors 

_A_(t)  *  [E^t),  E2(t)  .  .  .  En(t)]  (22) 

and  is  such  that_/^_(T)  =  _^V(0),  i.e.,  the  eigenvectors 
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are  periodic.  For  the  control  analysis  we  will  need  the 
eigenvectors  around  the  periodic  orbit  about  (ref.  7). 


Using  Equation  (10)  and  Equation 

( 21 ) ,  we  get 

A(t)|(t)  =  A(t)_A_(t)eJt_A_-1(0) 

(23) 

differentiating  (21)  leads  to 

A.eJtA._1(0)  +  A.^  _/Y-1<o> 

=  A(  t  )<§ 

(24) 

=  A_/\_  ( t )  e‘rV/\_-1  ( 0  ) 

(25) 

This  simplifies  to 

j±t)  -  A(t)_/\_(t)  -  _Ajt)J 

(26) 

This  is  the  equation  that,  when  integrated,  gives  the 
eigenvectors  around  the  periodic  orbit. 

The  Floquet  Theory  described  above  was  used  to 
investigate  the  stability  of  the  orbit  about  L3.  The  orbit 
was  found  to  be  unstable  since  one  of  the  system's  Poincare 
exponents  is  a  positive  real  root,  while  its  conjugate  is 
negative.  The  remaining  exponents  are  purely  imaginary. 
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4 


Table  I 


Poincare  Exponents  of  L3  Periodic  Orbit 

Mode  Poincare  Exponent 

Planar  0.0+0.77738i 

0 .0-0.77738i 


Planar 

Out  of  Plane 


0. 17893+0. Oi 
-0. 17893+0. Oi 

0 .0+0 . 84134i 
0.0-0.84134i 


In  Shelton's  orbit,  the  unstable  exponent  was 
designated  X3  and  is  the  root  that  his  controller  affected. 
Implementing  the  controller  yields 

=  A(T)ix(T)  +  B(t)u(t)  (27) 


where  ti(t)  is  the  control  vector  and  B(t)  is  the  control 
matrix.  Since  it  is  only  desired  for  the  controller  to 
affect  only  one  mode,  it  is  necessary  to  uncouple  the  modes. 
This  can  be  done  by  a  similarity  transformation 
(ref.  8:170-172).  Defining  the  uncoupled  variable  vector  as 
rj"(t),  we  can  obtain  rj”(t)  from 

fj’(t)  *  ^V~X(t)  4x(t)  (28) 

where  is  the  inverse  of  the  eigenvector  defined  in  (22). 

Differentiating  Equation  (28)  with  respect  to  time  yields 


& 


_At)ij  (t)  +  /^(tJrpt) 


(29) 


11 


W  } 

In  order  to  obtain  a  differential  equation  for  the  modal 
coordinates,  equate  Equations  (27  and  (28)  to  get 

*}  *  -2V  fA(t)^\_(t)  -_/^t>l  +  i|~1(t)B(t)u(t)  (30) 

From  (26)  we  had 

_A_=  a ( t >  _A_(t)  -  _/\_(t)j 
substituting  into  (30)  yields 
*—  —  —1  _ 

*|  ■  Ji|  +  .A.  (t)B(t)u(t)  (31) 

Since  the  system  is  diagonalized,  it  is  now  a 
relatively  simple  task  to  feed  back  the  single  unstable 
mode.  This  is  the  method  of  modal  control.  Since 
was  designed  as  the  unstable  exponent,  this  requires  that 
u  =  krj3.  This  yields 

"  J*\  +  ^V-1(t)B(t)Ki|3  (32) 
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IV.  Perturbation  Analysis 

Development  of  Perturbations 

In  the  restricted  four-body  problem,  i.e.  no 
inclination  or  eccentricity  of  the  moon  is  included,  the 
perturbation  from  the  reference  periodic  orbit  becomes 


where  X  are  the  actual  states  and  X  „  are  the  states 

p.  o. 

of  the  reference  periodic  orbit.  The  &x_  desig- 
nates  that  this  is  a  free  perturbation,  i.e.  no  forcing 
function  from  the  Moon's  eccentricity  or  inclination. 

Shelton  found  this  is  an  unstable  yet  controllable  system 
with  <&F  =  A(t)  £xp  (34) 

This  is  also  a  linearized  system,  Since  second  order  terms 
and  higher  were  truncated. 

In  the  unrestricted  four-body  problem  the  Moon's 
eccentricity  and  inclination  must  be  included.  Choosing  to 
define  the  perturbation 

£x  =  £xp  +  £xp  (35) 

where  &x_  is  defined  as  in  (33)  it  is  seen  that  &x 

F  p 

must  include  the  perturbation  due  to  the  forcing  function  from 
the  Moon's  eccentricity,  inclination  and  any  nonlinear  terms 
not  included  in  the  restricted  four-body  problem.  It  should 
be  noted  the  instabilities  in  the  restricted  four-body 
problem  arise  from  the  linearized  forcing  terms  due  to  the 
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inclusion  of  the  sun.  Therefore,  the  nonlinear  terms  in  the 
unrestricted  four-body  problem  are  the  nonlinear  part  of  the 
forcing  function  from  the  inclusion  of  the  sun  in  the  system 
model . 

Since  the  total  perturbation  is  derived  as 


ix  -  XR 


(36) 


where  XR  are  the  states  in  the  unrestricted  model, 
by  equating  (35)  and  (36)  it  is  seen  that 


=  X„  -  X  _  -  &x 

F  R  p.o.  p 

substituting  (37)  into  (34) 


-  A<t>  IXR  -  Xpo 


rearranging  (35)  yields 
£xp  -  £x  -  £xf 
substituting  (35)  gives 


(37) 


(38) 


(39) 


(40) 


Time  differentiating  (40)  and  substituting  (38)  produces 


<kc  =  XD  -  ^  o  -  A(t)  [XD  -  x  „  -  ixl  (41) 

p  R  p.O.  R  p.O.  p 

which  is  an  explicit  differential  equation  for  the  nonlinear 
perturbation . 
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Perturbation  Verification 

To  insure  that  the  states  and  perturbations  were 
being  calculated  correctly,  a  check  was  done  by  comparing 
the  states  derived  from  the  total  nonlinear  equations  of 
motion  with  those  derived  from  Equations  (4),  (37),  and 
(41).  This  was  done  in  the  following  manner: 

XD(t)  =  (t)  +  £x(t)  +  (42) 

R  p.O.  p  F 

But  from  Equations  (8),  (9),  and  (23) 

£xp(t)  *  _A_<t)  eJt  -1(0)  £xp(0)  (43) 

which  gives 

XR(t)  =  X  (t)  +  _A.(t)eJt  p  ( 0 ) 

+  ^Xp(t)  (44) 

using  Equations  (4),  (37),  and  (41)  for  _  (t),  &x„ 

p  •  o  •  r 

and  ^Xp(t),  respectively,  and  integrating  (4)  and  (41), 
the  state  vector  X  was  calculated.  This  state 

K 

vector  was  checked  against  the  state  vector  derived  from  the 
nonlinear  equations  of  motion,  (4),  and  was  found  to  agree 
extremely  well.  This  verified  the  states  and  perturbations 
were  being  calculated  correctly. 

With  Equations  (34)  and  (41)  the  controller  can  now 
be  implemented  in  such  a  way  as  to  affect  only  the  original 
instabilities  which  manifest  themselves  in  &xp. 
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and  by  ignoring  those  arising  in  &Xp,  which  if  they 
remain  small,  are  of  no  consequence  and  can  be  ignored.  In 
this  way  the  system  can  be  stabilized  in  the  "real"  dynamics 
of  the  complete  unrestricted  four-body  problem.  However,  to 
insure  non-excitement  of  the  stable  modes,  the  controller 
must  be  implemented  in  the  diagonalized,  uncoupled  modal 
coordinate  system. 
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V.  Development  of  the  Controller 


Control  Term  Derivation 

It  is  desired  to  control  only  the  "free"  part  of  the 
perturbation  ixF,  but  this  must  be  done  in  the  modal 
coordinate  system  to  insure  non-excitement  of  the  stable  modes. 

We  must  first  insure  that  the  modal  coordinates  can 
be  split  into  their  forced  and  free  components,  as  the 
perturbations  are  in  Equation  (35).  Equating  (35)  and  (28) 
produces 


*  A-1  <^xp  +  (46> 

defining 

-/V-1  <kp  -  *fp  (47) 

and  A.”1  <^F  =  *j"F 
Equation  (43)  becomes 

T  *  V  +  V  <48> 

and  it  is  seen  that  the  modal  coordinate  can  be  defined  in 
the  same  manner  as  Equation  (35). 
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Differentiating  the  free  part  of  Equation  (47) 

yields 

<^F  =  A.*fF  +  Xj\  F  (49) 

Equating  this  to  Equation  (27) 

=  -Mf  +  A/[f  (50) 

=  A(t)  ixF(t)  +  B(t)u(t) 
and  solving  for  rj’p 

rj'p  =  _A_“1(t)  [ A( t )  _A_(t)  -  _A_(T)]np  +  J/C1B(t)u(t)  (51) 

from  (26)  and  (32) 
is  obtained 

rj'p  *  jij"F  +  _A__1(t)  B(t)  Ku  (52) 

where  u  =  Krj  3  and  K  is  the  control  gain.  The 
B  matrix  defined  such  that  B3  =  [0,  1.,  0,  1.,  0,  0]T  (53) 

where  3  designates  the  third  column,  which  is  the  same  as 
used  in  previous  studies.  The  other  5  columns  are  zero. 

This  corresponds  to  two  fixed  thrusters  operating  at  a  45° 
angle  to  the  Earth-Moon  line. 
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In  a  similar  manner  the  equation  of  motion  for  the 


forced  perturbation  can  be  found  to  be 

Tp  “  -  Vo.  ‘  h(t) 

-  _/Y_1(t)  [A(t)_^_(t)  -  ifp  (54) 

It  is  necessary  to  implement  the  controller  in  the  real 
system  model,  since  a  controller  can  only  affect  the  "real 
world"  dynamics.  In  the  most  general  sense,  the  control 
vector  is  u  =  Krj'p  f  (55) 

substituting  for  rj"F  from  Equation  (28) 

u  -  K  JA_"1(t)^p(t)  (56) 

The  control  term  which  is  added  to  the  equations  of  motion 
( 4 )  becomes 


C.T.  *  Bu 

=  B  K  _A_“1(t)  £xp(t)  (57) 

using  the  definition  for  &x_  of  Equation  (37),  the 
control  term  becomes 

C.T.  =  B  K  _A__1(t)  [XR  -  Xp  Q  -  £xp]  (58) 

Now,  integrating  Equation  (41),  and  choosing  initial 
conditions  for  £xp»  the  control  term  can  be  calculated 
around  the  orbit  and  added  to  the  equations  of  motion  (4). 
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Calculation  of  Controller  Gain 


Equation  (31)  can  be  expanded  to 


11 

12 

13 
iA 


*1 2 
*1  3 
*1  4 

Is 

1« 


+  JC1  BU 


(59) 


Since  the  third  modal  coordinate  has  been  designated 
as  the  unstable  mode,  we  wish  to  feed  back  i^.  This 
resulted  in  choosing 


u  =  K  rj  3  ( 60  ) 

The  B  matrix  is  chosen  to  be  a  constant  matrix  and 
is  representative  of  the  physical  constraints  of  the 
problem.  Since  only  the  inertial  velocity  components  of  the 
satellite  can  be  controlled,  the  B  matrix  is  chosen  as 
defined  in  the  last  section.  This  corresponds  to  two 
reaction  jets  aligned  oppositely  and  45°  from  the  X-Y  axes. 
The  third  modal  coordinate  differential  equation,  with  the 
controller,  now  becomes 

rj3  *  U3  +  Ke(t)]  i|3  (61) 
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where 


€.  *  _A__1  ( t ) B  (62) 

Equation  (57)  can  be  solved  using  an  integrating  factor. 

€  is  available  as  a  Fourier  Series.  Separating  e3 

into  its  constant  and  sinusoidal  parts,  the  solution  to  (61) 

is 

r^(t)  ■  i|3(0)exp(\3  +  K63(,)t  exp  j^Ke3p(t)cxt  (63) 

Since  the  last  exponential  is  a  periodic  function  of  time, 
the  equation  for  the  new  Poincare'  exponent  for  the  third 
mode  simply  becomes 

\3  =  X3  +  K  <_A_~32  +  _A_~34 )  (64) 

using  the  value  for  the  old  Poincare  exponent,  \3,  and  the 
D.C.  part  for  ^V~32  and  _A_~34'  Equation  (64)  becomes 

\3  =  .1789305466  -  12 . 96306103 (K)  (65) 

Using  Equation  (65)  the  gain  K  could  be  chosen  to  place  the 
new  Poincare'  exponent,  in  the  stable  left  half  plane. 
Choosing  a  new  Poincare'  exponent  of  *  -.1789305466, 
the  gain  becomes  K  =  .027606218. 
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orbits,  with  no  moon's  eccentricity  or  inclination  included  and 
with  no  control  i.e.  K  =  0,  for  4  orbits.  As  can  be  seen  in 
Figures  4  and  5  r|^  and  ijj  exhibit  the  expected 
general  behavior.  From  Figure  6  it  is  seen  that 
grows  exponentially  as  expected.  However,  the  modes  grow 
until  they  exit  the  region  where  the  linearization  is 
valid.  Mode  four  verifies  this  behavior,  since  should 
decrease  exponentially,  which  it  does  initially  but  then 
oscillates  erratically  as  the  linear  regime  is  exceeded. 

Figures  8  through  11  show  the  modal  coordinates  versus  orbits, 
with  no  control,  but  with  eccentricity  and  inclination 
included.  It  is  seen  that  the  coordinates  almost  immedi¬ 
ately  are  within  the  nonlinear  region  (note  the  scale  differ¬ 
ence  from  Figures  4-7)  and  hence  the  system  is  totally 
unstable. 

To  verify  that  the  controller  was  working  correctly, 
the  Shelton  model  was  duplicated  and  the  controller  applied 
to  check  and  insure  the  model  could  be  controlled  as 
Shelton's  controller  had  done.  The  model  was  duplicated  by 
removing  the  nonlinear  terms  and  inclination  and  eccentri¬ 
city  effects  from  the  control  term.  This  results  in 
Equation  (58)  becoming 

C.T.  -  B  K  ^V"1  IXj,  -  Xp  0  ]  <66> 

which  indicates  that  the  controller  is  acting  only  on  the 
linear  perturbation  terms  in  the  restricted  four-body 
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Evaluation  of  the  Controller 

The  controller  was  designed  based  on  a  linearized 
set  of  equations  but  was  applied  to  the  full  nonlinear 
equations  of  motion.  Since  it  is  required  to  control  a 
satellite  with  all  effects  included,  the  controller  was 
included  in  the  full  nonlinear  equations  which  represent  the 
most  realistic  model  available. 

From  Chapter  IV  of  this  report,  without  the 
controller,  it  is  expected  that  and  rj2  will 
be  purely  oscillatory,  rj  3  will  exponentially  increase 
r|4  will  exponentially  decrease.  The  last  two  modes 
remain  zero  since  no  initial  condition  was  given  and  no 
forces  are  present  which  will  generate  these  components. 

This  can  be  seen  from  Equations  (14)  and  (31).  With  no 
control  (31)  becomes 


From  Table  I  we  can  see  that  the  modal  coordinates  will 
behave  as  described  above.  Figures  4  through  7  show  the 
modal  coordinates  n^  through  n4  versus  the  number  of 
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problem,  since  the  free  perturbation  term  is  defined  as  in 
Equation  (37).  Figures  12  through  15  show  the  results  of 
this  verification.  Figures  12  and  13  show  that  modes  one 
and  two  are  oscillatory  as  expected,  and  Figures  14  and  15 
reveal  that  modes  three  and  four  decrease  exponentially.  An 
exponential  curve  fit  was  done  on  mode  three,  the  controlled 
mode.  The  resulting  equation  was 

r|3  =  -.229e"*l61(t>  (67) 

The  Poincare  exponent  for  mode  three  is 

P.E.  =  -.179 

and  it  is  seen  that  the  two  exponents  agree  well  within  the 
accuracy  of  the  curve  fit. 
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VI.  Results,  Conclusions.  Recommendations 


Inclusion  of  Nonlinear  Perturbations 

The  next  step  in  the  control  analysis  was  to  include 

the  nonlinear  terms  from  the  restricted  four-body  problem  in 

the  model,  and  have  the  controller  act  on  not  only 

&x_,  but  the  as  well. 
r  p 

The  control  term  is  defined  in  Equation  (58)  as 

C.T.  -  B  K  Tl'1  (Xr  -  Xp  e  -  4ip] 

Figures  16  through  19  show  the  results  of  the  controller 
acting  in  this  model  with  £xp  included.  No  effects 
from  the  moon's  Keplerian  eccentricity  or  inclination 
have  been  included  in  this  model,  only  the  above  mentioned 
nonlinear  terms  from  the  sun  and  moon's  free  pertubations. 

As  can  be  seen  in  Figure  18,  modal  coordinate  three,  r^3, 
is  diverging.  When  rj  3  is  so  large  as  to  be  outside  the 
controllable  region,  all  the  modes  go  unstable,  as  is  indicated 
in  Figure  19  by  rj  ^  approaching  zero,  then  diverging 
after  approximately  three  orbits.  A  new  method  had  to  be 
devised  to  stabilize  not  only  the  free  perturbation,  &x_, 

r 

but  also  the  nonlinear  perturbations. 
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Encke's  Method 


The  method  used  in  designing  the  controller  to 
stablize  the  linear  and  nonlinear  perturbations  was  Encke's 
Method  (ref.  9:  390-396).  Shelton's  controller  had 
stablized  the  free  linear  perturbations  resulting  in  a  new 
reference  periodic  orbit.  Figure  20  shows  this  reference 
periodic  orbit  and  the  actual  orbit  the  satellite  would 
traverse.  The  difference  between  these  two  orbits  is  the 
perturbation  caused  by  the  nonlinear  effects.  In  Encke's 
Method,  when  the  perturbation  grows  larger  than  a  specified 
value,  rectification  occurs.  Rectification  means  that  a  new 
epoch  and  starting  point  is  chosen  which  coincides  with  the 
actual  orbit  of  the  satellite.  The  controller  is  activated 
at  this  point  to  zero  out  the  perturbation,  and  place  the 
satellite  back  on  the  desired  reference  orbit.  This  results 
in  higher  control  costs,  but  is  a  necessary  function  in 
order  to  stabilize  the  satellite.  Rectification  was 
achieved  by  zeroing  out  the  nonlinear  perturbation,  which 
results  in  the  controller  being  activated.  This  can  be  seen 
from  Equation  (58).  At  time  of  rectification,  the  nonlinear 
term,  &xp,  has  grown  to  such  a  value  that  it 
is  approaching  in  magnitude  the  free  perturbation,  ixp, 
which  is  defined  in  Equation  (37).  This  results  in  the 
control  term,  Equation  (58),  being  of  small  magnitude. 

By  setting  (ix^  to  zero,  the  controller  is  "turned  on" 

and  reacts  to  what  it  now  "thinks"  is  the  free  perturbation. 
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This,  in  effect,  "dumps”  the  nonlinear  perturbation  into  the 
free  perturbation.  This  "dumping”  is  perfectly  valid  because 
how  the  perturbation  is  split  up  arbitrarily  and  the  controller 
reacts  on  the  total  perturbation  when  rectification  occurs.  In 
order  to  determine  when  to  rectify,  an  average  value  of  the 
perturbation  states,  i.e. 

i  I  ,  I  !  ,  I 

(68) 


&x  * 

p 

+ 

*Px 

+ 

+ 

1  avy  1  4 

was  calculated.  When  the  value  exceeded  an  arbitrary  value, 
rectification  occured. 


Results  of  Encke's  Method 

Using  the  same  initial  conditions  as  were  used  when 
verifying  the  controller  would  zero  out  the  free  pertur¬ 
bation  only,  Encke's  Method  was  used  to  zero  out  both  the 
free  perturbation  and  the  nonlinear  perturbations.  Figures 
21  through  24  show  the  results  of  this  rectification  scheme. 
Figure  23  reveals  that  the  controller  zeros  out  the  initial 
transient  and  then  acts  as  a  dead  band  controller.  Modal 
coordinates  one  and  two  exhibit  the  expected  oscillatory 
behavior.  Coordinate  three,  the  unstable  mode  coordinate, 
demonstrates  stable  behavior.  The  periodic  spikes  are 
generated  when  rectification  occurs,  activating  the 
controller. 

Shelton  had  succesfully  stabilized  the  free 
perturbation  for  approximately  2  years  or  24  months. 

Figures  25  through  28  show  the  controller  acting  for  a 
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period  of  5  years  (60  orbits).  It  is  seen  that 
rectification  successfully  stablizes  the  free  and  nonlinear 
perturbations  for  this  period  of  time,  and  it  appears  would 
continue  to  do  so  for  much  longer  periods  of  time.  Figures 
25  and  26  also  show  the  sun's  effect.  This  effect  shows  up 
in  the  one  year  frequency  of  this  oscillatory  mode. 

The  moon's  eccentricity  and  inclination  were 
included  in  the  next  step  of  the  analysis.  This  was  to 
determine  if  Encke's  Method  and  the  redefinition  of  the 
perturbation  could  successfully  control  the  instabilities 
which  arise  from  these  lunar  forcing  functions.  Figures  29 
through  32  show  the  results  of  including  inclination  and 
eccentricity  of  the  moon.  It  is  seen  that  all  the  modes 
have  gone  unstable  within  six  orbits.  The  perturbations 
grow  rapidly  and  enter  the  nonlinear  region  almost 
immediately.  Since  the  controller  stablized  the  nonlinear 
terms  in  the  restricted  four-body  problem,  this  implies  that 
the  perturbations  arising  from  these  terms  do  not  grow  at  an 
uncontrollable  rate  or  are  merely  oscillatory.  However,  the 
perturbations  from  the  moon's  inclination  and  eccentricity 
do  grow  at  too  fast  a  rate  for  the  controller  to  handle. 

It  is  also  possible  that  one  or  both  of  the 
oscillatory  modes  becomes  unstable  with  the  inclusion  of  the 
Lunar  inclination  and  eccentricity.  Since  the  real  mode  is 
the  only  mode  controlled,  an  oscillatory  mode  going  unstable 
would  eventually  render  all  the  mode  unstable. 
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Control  Costs 


Without  the  Lunar  inclination  and  eccentricity 
included,  the  controller  successfully  stablized  the  satel¬ 
lite,  as  discussed  in  the  previous  section.  A  limitation 
was  put  on  the  magnitude  of  the  average  nonlinear 

perturbation,  £x  ,  defined  by  (68).  To  analyze 
pavg 

the  effect  of  varying  this  limitation  a  control  cost  or  an 
average  control  acceleration  must  be  computed. 


Since  control  is  an  acceleration  (A  velocity). 


-£‘ 


control  acceleration  I  dt 


where  Av  is  the  change  in  velocity.  An  average  acceleration 
is, 

_  1  f t 

A  *  —  I  (control  acceleration |  dt  (70) 

t  Jo 

where  A  is  the  average  acceleration. 

The  average  acceleration  was  calculated  for  several 
test  cases.  The  controller  must  first  zero  out  any 
transients  in  the  initial  values  of  the  states  and  will  then 


act  as  a  dead  band  controller,  allowing  the  satellite  to 

drift  until  the  maximum  value  of  £x  is  reached. 

pavg 

This  leads  to  the  average  acceleration  being  relatively 
large  initially  and  then  settling  down  to  a  constant  value 
after  a  period  of  time  has  elapsed. 


Figure  33  shows  the  average  acceleration  versus 

orbits,  and  exhibits  the  expected  behavior.  The  periodic 

spikes  are  caused  from  rectification  occuring,  and  the 

controller  is  activated.  The  average  acceleration 

approaches  a  fairly  constant  value  after  36  orbits. 

An  analysis  was  done  to  determine  the  effect  of 

varying  the  limit  of  £x  on  the  average  control 

pavg 

cost.  Runs  were  made  for  36  orbits  (3  years)  and  the  results 

can  be  seen  in  Table  II.  If  the  limit  on  &T  is  set 

pavg 

too  large,  the  control  acceleration  is  very  large,  indicating 

that  the  satellite  "drifts"  far  from  the  reference  orbit  and 

when  rectification  occurs,  the  controller  must  expend  a 

large  amount  of  energy  to  drive  the  perturbation  to  zero. 

As  the  limit  on  <f»x  is  tightened,  the  average 

pavg 

acceleration  decreases  as  expected  since  the  controller  does 
not  need  to  exceed  as  much  energy  to  drive  the  perturbation 
to  zero.  However,  the  number  of  rectifications  per  orbit, 
which  is  equivilant  to  the  number  of  times  the  controller  is 
activated,  increases  slightly,  as  expected.  Figure  34  shows 
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the  average  acceleration  versus  the  limit  on  &x  ,  on 

pavg 

a  log-log  scale.  This  plot  indicates  the  trends  discussed 
above. 


Table  II 


ixp 

av*  max 

A 

Rectification 
Per  Orbit 

1  x  10"2 

Large  (107) 

1  x  10"3 

.87  x  CIO”4) 

2.08 

1  x  10"4 

.86  x  CIO”5 ) 

1.86 

1  x  10“5 

.76  x  <10"6) 

1.94 

1  x  10"6 

.14  x  (10"6) 

2.11 

1  x  10"7 

.38  x  <1<T7) 

2.5 

Recommendations 

The  controller  demonstrated  its  ability  to  stabilize 
the  nonlinear  terms  in  the  perturbations  caused  from  the 
influences  of  the  earth,  moon,  and  sun.  When  the  moon's 
eccentricity  and  inclination  were  included,  the  controller 
could  no  longer  stabilize  the  system.  The  controller  was 
derived  based  on  modal  control,  and  only  one  mode,  the  real 
unstable  mode,  was  controlled.  It  is  possible  that  the 
oscillatory  mode  becomes  unstable  with  the  inclusion  of 
eccentricity  and  inclination.  This  would  result  in  the 
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total  sytem  becoming  unstable,  since  no  control  is  placed  on 
the  oscillatory  mode.  It  may,  therefore,  be  necessary  to 
place  some  form  of  control  e.g.  damping,  on  the  oscillatory 
mode.  Unfortunately,  we  are  not  free  to  choose  what  the 
controller  may  act  on,  since  the  controller  can  only  act  on 
the  momentum  components  of  the  perturbation.  A  limitation 
is  consequently  placed  on  how  the  modes  may  be  controlled. 
Therefore,  a  method  must  be  developed  that  will  allow 
placement  of  the  individual  Poincare  exponents,  within  the 
physical  restrictions  of  the  problem. 
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Appendix  A 

«  Angular  velocity  of  the  rotating  frame;  the 
moon's  mean  inertial  angular  velocity  vector 

m^  *  mass  of  earth 

m^  =  mass  of  moon 

m^  =  mass  of  sun 

G  *  Universal  Gravitational  Constant  (set  equal  to  1) 

m_  =  Scale  factor  locating  the  inertial  point  of 

-  system 

M 
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La grange  Points  in  Earth-Moon  System 
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Figure  2.  Wiescl  Orbit  About  L3 
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Fioure  3.  Geometry  of  the  Four-Body  Problem 
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Figure  4:  Modal  Coordinate  l,_No  Eccentricity  Or 
Inclination, (ox(0)«0,  4  Orbits) 
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Figure  10s  Modal  Coordinate^ ,  With  Eccentricity  and 
Inclination^  ox+5,  4  Orbits) 
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Figure  20s  Reference  Orbit  With  Rectification 
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Figure  28:  Modal  Coordinate  4,  Controller  With  Recti 
fication,  (Ox+tO),  €0  Orbits) 
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Figure  29:  Modal  Coordinate  1,  Eccentricity  and 
Inclination  Included 
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